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ABSTRACT 


A  derivation  Is  presented  of  the  equation  governing  the  pressure  In 
o  thin,  flat  film  of  ideal  gas  under  isothermal  conditions,  when  the 
surfaces  bounding  the  film  are  in  relative  normal  and  tangential 
motion.  When  tangential  motion  is  absent,  the  pressure  equation 
reduces  to  a  nonlinear  heat  equation  ,  which  admits  of  very  few 
closed-form  solutions.  Various  approximation  methods  are  discussed, 
and  two  problems  involving  small  periodic  variation  of  the  gap  be¬ 
tween  parallel  plates  are  solved  by  perturbation  method  for  a  film  in 
which  fluid  Inertia  is  negligible. 
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ISOTHERMAL  SQUEEZE  FILMS 


INTRODUCTION 

Because  of  its  wide  technological  application,  the  theory  of  fluid 
film  lubrication  between  surfaces  in  relative  lateral  motion  has  been 
extensively  studied  since  its  formulation  by  Osborne  Reynolds 
almost  eighty  years  ago.  The  recent  development  of  air -lubricated 
bearings  has  stimulated  the  extension  of  the  theory  to  the  case  of 
compressible  lubricating  films. 

Less  attention  has  been  given  to  the  pressure  generated  in  a 
fluid  film  by  relative  normal  motion  of  the  surfaces.  Most  of  the 
published  work  on  this  subject  has  been  confined  to  the  study  of 
incompressible  films  between  parallel  surfaces  and  in  journal 
bearings  QJ  .  Gas  squeeze  films,  as  they  are  called,  have  remained, 
for  the  most  part,  a  curiosity. 

The  purpose  of  this  paper  is  to  derive  and  apply  the  equation 
governing  the  pressure  in  a  thin,  flat  film  of  ideal  gas,  under  iso¬ 
thermal  conditions,  when  the  surfaces  bounding  the  film  are  in  relative 
normal  motion. 


In  a  recent  paper,  Elrod iZ/ derived  the  equation  governing  the 
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Bteady-Btate  pressure  in  a  journal  bearing  lubricated  by  an  incom¬ 
pressible  fluid  with  constant  viscosity.  Although  such  derivations 
are  usually  carried  out  with  the  help  of  several  ad  hoc  order-of- 
magnitude  assumptions,  Elrod  used  a  perturbation  approach,  with 
the  ratio  of  film  thickness  to  bearing  length  as  the  small  parameter. 
This  approach  not  only  allowed  him  to  derive  Reynold's  lubrication 
equation  in  a  more  convincing  manner,  but  also  enabled  him  to  retain 
the  terms  resulting  from  the  film  curvature.  Elrod  found  tliat  the 
Reynolds  equation,  with  the  first  order  correction  for  film  curvature 
included,  can  be  written 

in  which 

D  ss  Shaft  diameter, 
h  =  film  thickness, 

P  =  fluid  pressure, 

V  =  shaft  surface  velocity, 

X  =  distance  around  shaft  in  direction  of  rotation, 
z  =  distance  parallel  to  sliaft  axis, 

^  =  fluid  viscosity. 

Elrod's  result  led  a  reviewer  [3]  to  conclude  that  only  the  film 
curvature,  and  not  the  fluid  inertia,  has  an  influence  in  the 
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approximation  of  the  second  order.  However,  the  importance  of 
inertia  depends  upon  the  film  Reynolds  number,  which  does  not  appear 
in  Elrod's  result.  His  dimensionless  formulation  is  aimed  at  deter¬ 
mining  the  effect  of  film  geometry,  not  the  effect  of  film  dynamics. 

By  way  of  example,  his  normalized  pressure  'if  is  defined  by 

_ 

where  p  is  the  lubricant  density  and  hQ  is  a  typical  value  of  the  film, 
thickness.  Since  the  pressure  in  a  journal  bearing  film  varies 
according  to  h^"^,  this  formulation  is  quite  adequate  for  Elrod’s 
purposes.  However,  to  represent  properly  the  dependence  of 
pressure  on  the  parameters  of  the  fluid,  as  predicted  by  lubrication 
theory,  it  is  correct  to  use 

-rr'-  ^  - 

which  differs  from  Elrod's 'IT'precisely  by  the  factor  — 

a  Reynolds  number.  That  inertia  effects  are  negligible  is,  then,  a 
tacit  assumption  made  by  Elrod,  rather  than  a  conclusion  to  be 
drawn  from  his  results. 

In  Sec.  I  of  the  present  paper,  we  use  an  approach  quite  similar 
to  Elrod's  in  order  to  derive  the  equation  governing  the  pressure  in 
a  thin  film  when  both  lateral  and  relative  normal  motions  of  the 
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surfaces  are  present.  We  consider  only  flat  films:  i.  e.  ,  although 
the  surfaces  bounding  the  film  may  be  curved,  it  is  possible  to  choose 
a  plane  such  that  the  distance  from  any  point  in  the  film  to  the  plane 
is  small  compared  with  the  lateral  dimensions  of  the  film.  The  analy¬ 
sis  does,  however,  apply  to  Journal  bearings,  subject  to  curvature 
corrections  of  the  type  derived  by  Elrod.  We  assume  that  the  density 
of  the  lubricant  is  proportional  to  its  pressure,  a  condition  which  is 
met  when  the  lubricant  is  an  Ideal  gas  and  isothermal  conditions 
obtain.  Although  fluid  inertia  is  usually  negligible  in  lubricating 
films,  we  do  not  assume  this  a  priori.  Rather,  we  derive  criteria 
under  which  inertia  can  validly  be  neglected,  and  indicate  a  method 
of  procedure  which  can  be  followed  when  these  criteria  are  not  met. 

In  Sec.  II,  we  introduce  the  restriction  that  lateral  motion  be 
absent,  and  thereafter  consider  only  pure  squeeze  films. 

A  few  exact  solutions  to  the  squeeze -film  equation  are  presented 
in  Sec.  III.  In  Sec.  IV,  we  consider  limiting  forms  of  the  squeeze- 
film  equation  at  high  and  low  frequency  of  the  squeeze  motion.  We 
find,  in  particular,  that  at  low  frequencies,  the  pressure  in  an  iso¬ 
thermal  squeeze  film  is  governed,  to  first  approximation,  by  the 
incompressible  squeeze -film  equation. 

In  Sec.  V,  we  calculate  the  pressure  field  generated  by  small 
periodic  variation  of  the  gap  between  infinitely  long  parallel  plates, 
and  in  Sec.  VI,  we  treat  the  equivalent  axisymmetric  problem. 
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The  methods  of  analysis  of  Secs.  V  and  VI  can  be  readily 
adapted  to  other  classes  of  problems  in  squeeze-film  theory.  The 
treatment  of  films  bounded  by  surfaces  other  than  parallel  flat  plates 
is  more  complicated,  but  not  conceptually  more  difficult.  Also,  cases 
involving  both  lateral  motion  and  squeeze  motion  can  sometimes  be 
solved  by  analogous  methods.  Current  research  in  the  IBM  General 
Products  Division  San  Jose  laboratories  includes  problems  in  each  of 
these  categories. 

W.  A.  Michael,  of  this  laboratory,  has  developed  an  alternate 
approach  to  the  solution  of  problems  in  gas  squeeze-film  theory.  This 
approach,  in  which  the  pressure  is  written  as  a  Fourier  series  in  the 
spatial  variables  with  time -dependent  coefficients,  has  been  used  by 
Michael  to  study  rectangular  squeeze  films.  His  report,  currently  in 
preparation,  includes  solutions  both  for  periodic  motion  and  for  the 
transient  response  to  a  sudden  variation  of  the  film  thickness. 
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I.  THE  MATHEMATICAL,  FOUNDATIONS  OF  ISOTHERMAL 

GAS -FILM  THEORY 


Assume  that  a  thin,  continuous  film  of  ideal  gas  is  contained 
between  the  surfaces 


H  (x.^  i)  , 

(1.1) 

-t) ^ 

(1.2) 

where  xj,  X2,  X3  are  right-handed  Cartesian  coordinates  and  t  is 
time.  The  film  thickness  h,  defined  by 

is  positive  for  all  values  of  xj^,  X2.  t;  the  surfaces  move 
relative  to  the  ambient  gas  with  velocity  components 

The  surfaces  bounding  the  film  may  be  either  rigid  or  flexible, 
but  are  assumed  continuous.  At  each  point  of  each  surface, 
three  components  of  velocity  provide  one  degree  of  freedom  too  many, 
and  must  therefore  be  related  through  a  kinenmtic  constraint.  With 
the  convention  that  Latin  indices  extend  over  the  values  1,  2,  and 
Greek  indices  over  1,  2,  3,  the  kinematic  constraints  are 


\l'-.  . 


repeated  indices  denoting  summation. 


The  motion  of  the  gas  in  the  film  is  governed  by  the  equations 
of  viscous  hydrodynamics.  The  viscosity  coefficients ^  and  X 
are  assumed  constant,  since,  within  the  range  of  interest  of  gas 
lubrication  theory,  their  variation  is  slight.  As  a  consequence,  we 
can  write  the  Navier-Stokes  equation 
Di/g 


P 


Dt 


<>2^1 


(1.6) 


in  which  yo  and  denote,  respectively,  the  density  and  pressure 
of  the  gas,  and  D/Dt denotes  material  differentiation; 

0/of  =  Vdt  +  • 

The  gas  also  obeys  the  continuity  equation 

where  the  dilation  A  is  defined  by 


(1.7) 


(1.8) 


A= -L 

P  Dt  • 

Thus,  Eq.  (1. 6)  can  also  be  written 


(1.9) 


Dt 


(1.  10) 


Since  the  film  thickness  is  small  compared  with  the  bearing 
breadth  B,  the  hydrodynamic  equations  can  be  simplified  by  using 


(1.11) 


I 


as  a  perturbation  parameter,  where  hQ  is  a  typical  value  of  the  film 
thickness.  Using  the  technique  employed  by  Elrod^ZJ  to  study  curved 
incompressible  films,  we  introduce  a  dimensionless  coordinate  system 
that  "stretches"  the  coordinate  normal  to  the  film: 

Xi  =  .  (1. 12) 

2  =•  .  (1.1.3) 

The  ^=1,2  components  of  Eq.  (1.  10)  then  become 


M  — i-  =  6 


+  6 


Dt 


'  (1.14) 


and  the  K  =  3  component  becomes 


We  now  introduce  dimensionless  variables  in  such  a  way  that  the 
magnitude  of  each  term  in  the  equations  of  motion  is  represented  by 
its  coefficient,  bearing  in  mind  that  we  intend  to  apply  the  equations 
to  a  lubricating  film. 

For  the  time  scale,  we  choose  the  reciprocal  of  a  typical  fre¬ 
quency  of  the  squeeze  component  of  surface  motion.  Thus,  we 


assume 


V  =  lloiXjW,  \/'=  ho^^Vl' 


(1.16) 
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where  W  and  W  are  dimensionless  velocities  of  order  unity,  and 
introduce  a  dimensionless  time  T,  defined  by 

T=  ojt.  (1.  17) 

To  be  consistent  with  Eq.  {1.  16),  we  let 

\f^  =  hoUJ^.  (1*18) 

The  scale  of  the  lateral  velocity  components  is  not  necessarily 
related  to  the  scale  of  the  squeeze  component.  While  it  is  true  that 
the  squeeze  motion  forces  gas  outward  or  sucks  it  inward,  at  a  char¬ 
acteristic  velocity  t<J  B,  there  is  also  a  contribution  to  the  lateral 
velocity  arising  from  the  lateral  motion  of  the  bearing  surfaces.  It 
is  always  possible  to  choose  an  instantaneous  orientation  of  the  coordi¬ 
nate  system  such  that  the  components  of  surface  motion  in  the  xi  and 
directions  are  of  the  same  order.  We  shall  assume  that  the  time  variation 


and  expect  the  to  be  of  order  unity. 
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In  terms  of  our  dimensionless  quantities,  the  constraints  (1.4) 
and  (1.5)  become,  respectively. 


Ul> 


iT  s 


U>VJ  ‘ 


OJ 


B 


j  H 

(1.21) 

(1.22) 

in  which 


^‘/^o  ■  (1-23) 

The  factor  V/^B  implicit  in  Eqs.  (1.21)  and  ( 1- 22)  provides  a  measure 
of  the  relative  magnitude  of  the  two  reciprocal  times  characteristic 
of  the  bearing  kinematics:  V/B  represents  a  shear  rate  character¬ 
istic  of  the  lateral  motion;  U)  ,  as  defined  above,  is  a  typical  frequency 
of  the  squeeze  motion. 

We  now  turn  our  attention  to  the  definition  of  a  dimensionless 
pressure.  The  theory  of  lubrication  for  incompressible  films,  which 
provides  at  least  a  limiting  description  of  gas  films,  shows  that,  in 
the  absence  of  squeeze  motion,  the  pressure  in  the  interior  of  a  film 
varies  according  to^^^V/fi^^  ^For  pure  squeeze  films  (lateral  motion 
absent),  the  variation  is  as^^Ctlt^.  Therefore,  when  lateral  and 
squeeze  motions  are  both  present,  we  may  find  useful  the  normaliza¬ 
tion 

ya/  =yu.  (ut  V-  ]/Iq j  ir^ 


(1.24) 
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steep  gradients  of  pressure  obtain,  and^j_i  may  be  quite  large. 

In  the  interior  of  the  film,  however,  they  will  normally  be  of  order 
unity.  An  obvious  exception  occurs  when  the  bearing  undergoes 
lateral  vibration  at  a  frequency  large  compared  with^.  In  a  step 
bearing,  *  be  large  In  the  neighborhood  of  the  step. 

Both^g  and^L  are  likely  not  to  be  of  order  unity  in  the  vicinity  of 
a  source,  unless  we  consider  only  the  completely  viscous  film  in 
which  fluid  inertia  is  negligible. 

We  now  Introduce  modified  Reynolds  numbers  corresponding, 
respectively,  to  the  squeeze  motion  and  to  the  lateral  motion: 

.  (1.31) 

In  terms  of  the  dimensionless  quantities  introduced  above,  the 
equationsof  motion  (1.  14)  for  the  lateral  velocity  components  become 

^  <j  S!s  _ 

.  0-i-  V/'^ B)  ^ 

(1.32) 

Equation  (1.  15),  for  the  normal  velocity  component,  becomes 

*  Although  we  have  been  considering  only  films  with  continuous 
bounding  surfaces,  piecewise  continous  surfaces  can  be  treated 


in  an  obvious  manner. 
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^  -  Tit  v/^s;  ^ 


(Rsi-f?i,)R«.iJw/jXi J-  ji*  ’’’  Yu-  v/u. 


(1.33) 


Equation  (1. 33)  implies  that,  with  neglect  only  of  terms  of  the 
second  degree  or  higher  in  the  pressure  is  constant  across  the 
This  conclusion,  of  course,  does  not  necessarily  apply  near 
the  periphery  of  the  bearing  or  in  other  regions  where  one  or  both  of 
the  dimensionless  dilational  stresses  become  large  (of  order 

^  *?  / — s. 

c  ).  Moreover,  either  PRs  or  PR^,  could  be  of  order  t  :  In  this 
unlikely  case  also,  the  pressure  may  vary  significantly  across  the 

In  most  cases  of  interest,  however,  it  is  correct  to  infer  from 
Eq.  (1.  33)  that  vanishes  throughout  the  interior  of  the  film, 

and  we  shall  proceed  on  the  assumption  that  this  is  the  case.  Con¬ 
sistent  with  this  assumption  is  the  reduction  of  Eq.  (1. 32)  to 

In  most  bearing  applications,  the  Reynolds  numbers  are  negli¬ 


gibly  small,  but  for  the  present  we  need  not  assume  this  is  so.  In 
view  of  Eqs.  (1.19),  (1. 20)  and  (1.  23),  the  dimensionless  velocity 
components  satisfy  the  boundary  conditions 
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UJ  -  = 

^  (i+  ojS/\/)  ' 


{/■/■  uBI\/)  ' 


(1-35) 


Since  'JJ',  and  consequently  P,  are  assumed  constant  across  the 
film,  Eq.  (1.34)  and  the  boundary  conditions  (1.35)  taken  together  are 
equivalent  to  the  integro -differential  equation 

27  i- 


H 

in  which 


(1.36) 


H  =  --  . 


(1.37) 

Equation  (1.36)  represents  a  preliminary  result  to  which  we  shall 
return  presently.  For  the  moment,  however,  let  us  consider  again 
the  continuity  equation  (1. 8).  In  terras  of  the  original  (dimensional) 


variables,  this  equation  can  be  written 
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<^p/iit  +■  'O,  (1-38) 


Since  we  are  neglecting  the  transverse  pressure  variation,  and  con¬ 
sequently  the  transverse  density  variation,  integration  of  Eq.  (1.  38) 
across  the  film  yields 

‘'^*1  +  f  ( Vi'-  %)-  0  ■  ».  39) 


However, 


so  that 


(1.41) 


In  view  of  the  kinematic  conditions  (1.4)  and  (1.5),  Eq.  (1.41) 


becomes 


Cpk)  _  _ 


(1.42) 


If  we  again  neglect  the  density  variation  across  the  film,  Eq.  (1.42) 
becomes,  in  terms  of  the  dimensionless  variables, 

_ ! _  j _ ^ _  ^TT  j  ~  ^?(1.43) 


(H  v/w  a; 
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Integrating  Eq.  (1.  36)  across  the  film  yields 

rF' 


_  _di 

J  Zif-hcuB/V)  /Z 


6W 


(1.44) 


-H')  [( ^ -h -t 


Substituting  Eq.  (1.44)  into  Eq.  (1.43)  yields 


— ^ - ^  Irr  H ( Ui  4- Uj:) 

(t-h  Vfcu6)  ojy  (!■¥  ooB/  v)  L 


(1.45) 


v/here 


H 

If.  =  ^  TT  P J  {2-  J  ( 1 . 46 ) 


r  1 

i,=  1,4  4^^  j  . 

^  ^  yV  “*  (1-47) 


We  note  that,  if  H,  Uf  and  are  specified  functions  of  Xj  and 
T,  the  dimensionless  pressure  is  the  only  dependent  variable 
appearing  on  the  left  side  of  Eq.  (1.45).  Consequently,  if  conditions 
are  such  that  the  right  side  can  be  neglected,  we  obtain  a  single 
partial  differential  equation  for  the  pressure  in  an  isothermal  gas 


film:  a  generalized  Reynolds  equation. 

The  right  side  of  Eq.  (1.45)  is  comprised  of  terms  arising  from 
the  inertia  terms  in  the  Navler -Stokes  equation.  The  inertia  factors, 
defined  by  Eqs.  (1.46)  and  (1.47),  involve  the  dependent  variables  TT> 
and  w,  and  hence  cannot  be  specified  a  priori.  However,  the 
various  dimensionless  quantities  have  been  defined  in  such  a  way  that, 
under  normal  circumstances,  Ig  and  are  of  unit  order.  Thus,  the 
significance  of  the  right  side  of  Eq.  (1. 45)  is  measured  solely  by  the 
modified  Reynolds  numbers  Rg  and  R^.  If  they  are  very  small  com¬ 
pared  with  unity,  Eq.  (1.45)  reduces  to 


describing  the  dynamics  of  the  bearing  system. 


Even  if  fluid  inertia  cannot  be  completely  neglected,  it  may  be 
possible  to  account  for  its  effect,  at  least  approximately,  if  the  modi¬ 
fied  Reynolds  numbers  Rg  and  R^  are  not  too  large.  Estimates  of 
the  inertia  factors  Ig  and  Ij^,  based  on  some  reasonable  assumption 
for  the  velocity  profile,  can  and  have  been  usedj^"^.  Equation  (1.47) 
thus  becomes  a  partial  differential  equation  with  the  dimensionless 
pressure  as  the  only  dependent  variable,  instead  of  an  integro- 
differential  equation  in  four  dependent  variables. 


II.  THE  SQUEEZE-FILM  EQUATION 

In  deriving  the  Reynolds  equation  (1.49),  we  found  it  convenient 
to  introduce  the  dimensionless  pressure  TT',  defined  by  Eq.  (1.24), 
since  this  quantity  is  normally  of  order  unity  in  the  interior  of  the 
gas  film.  In  this  section,  however,  we  are  concerned  with  the  Impli¬ 
cations  of  Eq.  (1. 49)  rather  than  with  its  derivation.  Since  the  load- 
bearing  ability  of  a  gas  film  is  measured  by  the  gage  pressure  in  its 
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interior,  it  is  useful  to  normaUze  the  pressure  with  respect  to  the 
ambient  pressure.  The  Reynolds  equation  then  becomes 


=  [pH(U.  +  Uj]+<r4^,  ,a.., 

Where  the  dimensionless  parameters /\  (the  bearing  number)  and 
(the  squeeze  number)  are  defined  by 

,  ,2.2, 

^  ^  ^  (2.  3) 

If  Uii  Uj[,  and  H  are  specified  functions  of  Xj  and  T,  Eq.  (2.  1) 
is  a  quasilinear  parabolic  differential  equation  in  one  dependent  and 
three  independent  variables. 


Spatial  symmetries  can  sometimes  be  used  to  eUminate  one  of 
the  Independent  variables.  The  simplest  example  of  this  procedure 
occurs  when  the  gas  film  is  infinitely  long  in  the  X^  direction  and  H 
is  independent  of  X^.  Omitting  the  subscripts  from  U^,  and  X^, 
we  reduce  Eq.  (2. 1)  to 


It  is  sometimes  possible  to  take  advantage  of  a  spatial  symmetry 
by  introducing  a  curvilinear  coordinate  system  Yt  Equation  (2.  1)  then 
becomes 


=  /I  [PH  ((iV  PH  I,  (U^+U^'J 


4-or  . 


where  and  denote  the  contravariant  components  of  the  lateral 
surface  motion,  ~  denote  the  components  of  the  contravariant  metric, 

n  ^ 

and  are  the  Cristoffel  symbols  of  the  second  kind. 

The  most  widely  used  curvilinear  coordinate  system  is  the  polar 


system 


&=aA.e:ia^  (  Xii/X,)^ 


(2.6) 


for  which  Eq.  (2.  5)  reduces  to 


=  Afe  [f?pn  ( u«t  [ph(  u;)  ]J  + 


(2.  7) 


in  which  Uj^,  U^,  denote  the  physical  components  of  surface 

motion.  In  the  important  case  of  axially  symmetric  motion,  we  obtain 


4  ( Rrf^P If  ) -  a|^  [fi  PH  c  +.  +  f?»- 


a  (PH)  , 


■pp“  •*  (2.8) 


when  the  bearing  surfaces  are  rigid,  the /\  -term  vanishes 


The  character  of  solutions  to  Eq.  (2,  1),  and  its  equivalents  (2.  5), 
(2.  7),  is  determined  by  the  magnitudes  of  and<r^.  In  the  remain¬ 
der  of  this  paper,  we  consider  the  special  case  of  zero  bearing  number 
Thus,  we  focus  our  attention  on  the  squeeze-film  equation 


_  aP  .  a(PH) 


(2.9) 


which  governs  the  film  pressure  in  the  absence  of  lateral  surface 
motion,  and  upon  its  axisymmetrlc  equivalent 


i(PH) 


(2.  10) 


For  infinitely  long  films,  Eq.  (2.  9)  reduces  to  the  one -dimensional 


squeeze-film  equation 


III.  EXACT  SOLUTION  TO  THE  SQUEEZE-FILM  EQUATION 


Since  the  squeeze-film  equation  (2.4)  is,  in  general,  a  nonlinear 
parabolic  equation  with  variable  coefficients,  present  day  analytical 
methods  can  be  expected  to  provide  exact  solutions  only  in  the  most 
degenerate  cases.  Fortunately,  one  of  these  cases  is  of  engineering 
interest;  Externally  pressurized  films,  under  steady  conditions,  are 
governed  by  Eq.  (2.  9)  with  the  squeeze  numberO —  set  equal  to  zero: 


(3.  1) 


On  source-free  segments  of  the  bearing  periphery,  the  pressure  is 
ambient;  at  supply  holes,  it  equals  the  inlet  pressure. 

Equation  (3.  1)  can  be  rewritten 


_^L 


=  0 


(3.  2) 


Since  Eq.  (3.  2)  is  linear  in  —  it  is  obtained  from  the  equation  for 
incompressible  externally  pressurized  films  simply  by  replacing  P 
with  P  —it  can  often  be  solved  by  straightforward  methods.  In 
particular,  if  the  externally  pressurized  bearing  consists  of  parallel 
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flat  plates,  so  that  H  is  constant,  Eq.  (3.  2)  reduces  to 


and  the  methods  of  potential  theory  become  available.  Because  of 
the  widespread  use  of  externally  pressurized  gas  bearings,  many 

theoretical  Investigations  of  their  properties  have  appeared  in  the 
literature. 

We  now  return  to  Eq.  {2.  9)  with  a  nonvanishing  right  side.  This 

iSi  in  essence,  a  nonlinear  heat  equation,  the  analogy  with  heat  flow 

2 

being  more  evident  if  we  think  of  P  as  the  dependent  variable.  This 
suggests  that  useful  results  might  be  found  in  Crank's  extensive  treatise 
On  page  162,  Cfank  describes  an  analysis  due  to  Wagner  which 
leads  us  to  an  exact  solution  to  Eq.  (2.  9).  Wagner  considered  diffusion 
of  a  solute  into  a  semi -infinite  medium,  with  the  diffusivity  propor¬ 
tional  to  the  concentration,  the  surface  concentration  held  constant, 
and  the  Initial  concentration  equal  to  zero. 

The  squeeze -film  analog  of  Wagner's  problem  is  the  determination 
of  the  pressure  field  between  parallel  plates,  originally  in  intimate 
contact,  then  suddenly  pulled  apart.  Actually,  we  can  generalize  the 
solution  to  include  an  arbitrary  step -function  change  in  the  gap  between 
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semi -infinite  parallel  plates.  Strictly  speaking,  this  problem  is  out¬ 
side  the  scope  of  squeeze-film  theory  as  reported  in  this  paper 
because,  at  the  initial  instant,  dllational  stresses  and  inertia  are 
surely  important,  and  the  caloric  behavior  of  tiie  film  is  probably 
nearer  adiabatic  than  isothermal  (if.  Indeed,  these  thermostatic 
concepts  have  any  meaning  at  all).  It  can  be  hoped,  however,  that 
squeeze-film  theory  provides  an  accurate  picture  once  the  first  few 
Instants  of  time  have  passed. 

We  assume  that  the  projection  of  the  gas  film  on  the  XJ-X2  plane 
occupies  the  half -plane  X  ^0.  Since  no  bearing  breadth  can  be  defined, 
we  normalize  all  spatial  variables  with  respect  to  the  final  gap  h^. 

We  suppose  that,  at  time  zero,  the  gap  is  suddenly  changed  from6)i.hp 
to  hg,  then  held  at  the  new  value.  Thus, 

L,i^o 

.  (3.4) 

U, 

We  use  the  lower  case  t,  denoting  actual  time:  since  there  is  no 
characteristic  frequency,  the  time  variable  cannot  be  normalized. 

For  positive  values  of  t,  the  normalized  pressure  P  is  governed  by 


where  o— *  is  a  characteristic  time,  defined  by 

Corresponding  to  the  discontinuous  change  in  gap  at  t  =  0  there 
is  a  discontinuous  change  in  the  gas  pressure.  With  the  isothermal 
conditions  postulated,  Boyle's  law  requires  that 

0)~  (3_7j 

Joyyv 
X-P  -0 


Since  the  pressure  is  ambient  at  the  bearing  edge^ 

^i0^i)^\  (3.9) 

Equation  (3.  5),  subject  to  the  initial  condition  (3.  7)  and  to  the 
boundary  conditions  (3.  8),  (3.  9),  admits  of  the  self -similar  solution 

(3.10) 

where  the  function  f  (y)  is  determined  by  the  ordinary  differential 


equation 
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subject  to  the  boundary  conditions 

~  4  Jorrv  /  (^)  ~  ec  ^ 

%  (3.12) 

Numerical  integration  of  Eq.  (3.  11)  subject  to  the  boundary  con¬ 
ditions  (3.  12)  yields  the  family  of  curves  illustrated  in  Fig.  3, 1,  As 
<k  becomes  arbitrarily  small,  the  solution  tends  uniformly  to  the 
curve  for  o4  =0  (plates  initially  in  contact).  This  limiting  solution 
joins  the  y-axis  at  about  y  =  0.81.  iUustrating  the  curious  fact, 
observed  by  Wagner  |67  .  that  Eq.  (3.5)  -  although  parabolic  - 
admits  of  a  solution  with  a  well-defined  wave  front.  Needless  to  say, 

the  basic  assumptions  of  squeeze -film  theory  cannot  hold  true  for  such 
an  extreme  case. 
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IV.  LIMITING  FORMS  OF  THE  SQUEEZE-FILM  EQUATION 

We  began  the  previous  section  by  considering  externally  pres - 
surized  films  at  zero  squeeze  number.  Nontrivial  solutions  result 
only  because  of  the  boundary  conditions;  The  pressure  at  supply  holes 
differs  from  the  ambient  pressure  at  source-free  segments  of  the 
boundary.  Were  this  not  so,  the  pressure  would  be  ambient  every¬ 
where  in  the  film,  as  is  physically  evident. 

When  the  squeeze  number  is  finite  but  small,  we  expect  that  the 
pressure  In  a  self-acting  squeeze  film  (no  external  pressurization) 
will  not  differ  much  from  ambient:  Small  squeeze  number  corresponds 
to  low  frequency,  so  that  the  bearing  has  time  to  "leak.  "  Thus,  with 
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iX.  Jj 

,(n-0^j 

Tf 


C>fX,3,... 


(4.4) 


Since  the  pressure  is  ambient  at  the  bearing  periphery,  all  the  p(*^) 


must  vanish  there. 


It  is  of  interest  to  let 


/+<rP“^=  /r. 


so  that,  with  Eq.  (4.  1), 


P»  ][’+  0  (<^*)  .  (4 

The  boundary  conditions  require  that  77"  be  unity  on  the  bearing 


periphery,  and  Eq,  (4.3)  implies  that 


j. 


This  is  the  incompressible  squeeze-film  equation,  for  which  many 
solutions  are  available  Qj  ,  Thus,  with  neglect  only  of  terms  of  the 
second  degree  or  higher  in  the  squeeze  number,  an  Isothermal  gas 
film  can  be  considered  incompressible.  Physically,  at  low  squeeze 
numbers  the  gas  leaks  out  before  it  is  significantly  compressed  or 
rarefied. 

At  very  large  squeeze  number,  on  the  other  hand,  the  gas  exhibits 
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almost  no  tendency  to  leak.  When  the  frequency  of  the  squeeze  motion 
is  sufficiently  high,  the  escape  of  the  gas  is  effectively  blocked  by  its 
own  viscosity:  The  bearing  tends  to  behave  like  a  bellows.  We  use 
the  perturbation  parameter  — ■  and  note  that,  as  the  squeeze  number 
becomes  indefinitely  large,  Eq.  (2.  9)  tends  to  the  limiting  form 

(4.8) 

There  is  obviously  something  wrong  with  this  equation:  It  predicts 
that  the  product  PH  is  a  function  of  the  spatial  variables  only,  which 
leaves  no  way  to  satisfy  the  boundary  conditions  on  the  bearing  peri¬ 
phery.  For  example,  if 

/  ^ 

so  that 

,  (4.10) 

Eq,  (4.8)  predicts  that  the  pressure  satisfies  Boyle's  law: 

P(X;,Tj-  i/|(T)  (J^O)  (4.H) 

throughout  the  Interior  of  the  film.  Thus,  the  pressure  must  drop 

(or  jump)  discontinuously  to  ambient  at  the  bearing  periphery.  Since 

-1 

Eq.  (4.8)  is  derived  from  Eq.  (2.9)  simply  by  setting  =  0,  so  that 
the  bearing  is  precisely  a  bellows,  it  might  be  hoped  that  the  next 


order  of  approximation  predicts  a  continuous  pressure  leakage  at  the 
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periphery,  but  this  is  not  the  case.  If  we  set 

(4.12) 

substitute  Eqs.  (4.  9)  and  (4.  12)  into  Eq.  (2.  9),  and  neglect  terms  of 

-1 

the  second  degree  and  Mgher  ing-'  ,  we  obtain 

^L'^'^)/^T=0  (4.13) 

with  the  same  difficulties  as  before. 

-1 

The  problem,  of  course,  is  that  perturbation  oner*  is  a  singular 
perturbation:  Eliminating  the  term  involvingCT"  from  Eq.  (2.  9) 
reduces  its  order.  It  would  seem,  then,  that  the  pressure  does  indeed 
obey  Boyle's  law,  except  that  in  narrow  boundary  layers  near  the  peri¬ 
phery,  it  changes  steeply  but  continuously  to  ambient.  Even  this  con¬ 
cept  must  be  treated  cautiously,  however,  for  the  bearing  leakage  is 
a  continuing  process:  The  boundary  layers  tend  to  diffuse  away,  as 
illustrated  in  Sec.  IV  by  the  exact  solution  for  the  semi- infinite  film. 

Mathematically,  the  problem  is  b^'ofold.  First,  passing  from 
Eq.  (2.  9)  to  Eq.  (4.  8)  involves,  as  far  as  spatial  derivatives  are  con¬ 
cerned,  a  drop  in  order  of  2,  not  1.  Also,  in  (Xj,  Xg.  T)  space,  the 
boundary  layers  run  parallel  to  the  characteristics  of  the  reduced 
'equation.  Both  difficulties  tend  to  obstruct  a  singular  perturbation 
approach.  The  problem  invites  attention. 


-  31  - 


V.  SMALL.  PERIODIC  VARIATION  OF  THE  GAP  BETWEEN 
INFINITELY  LONG  PARALLEL  PLATES 


Consider  an  infinitely  long  bearing  of  uniform  width  B.  We  can 
choose  a  Cartesian  coordinate  system  so  that  the  projection  of  this 
bearing  on  the  X3  =  0  plane  is 


I  • 

I  —  00  ^  ^  oa  J 

The  one-dimensional  squeeze-film  equation  (2.  11)  is  then  applicable. 

Assume  that  the  surfaces  of  the  bearing  are  parallel  plates  and 
that  the  gap  between  them  varies  according  to 


h  =  ho  (/  + cot)^ 


where  the  constant  S  is  small  compared  with  unity.  We  then  expect 
that  the  deviation  of  pressure  from  ambient  will  be  of  order  ^ 

Thus, 


We  now  substitute  Eqs.  (5.2)  and  (5.3)  into  Eq.  (2.  11)  and  neglect 
terms  of  the  second  degree  or  higher  in  .  In  terms  of  the  normalized 
variables,  we  obtain 


(5.4) 


-  32  - 


Since  the  pressure  is  ambient  at  the  bearing  edges, 


7r(±/jl,T)=  0  . 


We  seek  a  steady-state  solution  to  Eq.  (5.  4)  subject  to  Eq.  (5.  5) 


by  assuming 


rrUj)=  Tr-U]&^r-h  'rQ(x)^T. 


Substituting  into  Eq.  (5.4),  we  obtain 


(<r+<r7I^  iaXrv7~_  (5  7) 

If  this  equation  is  to  hold  for  all  values  of  T,  both  sides  must 
vanish  identically.  Similarly,  the  boundary  conditions  expressed  by 
Eq.  (5.  5)  are  satisfied  if  and  only  if  both  and  TJ^vanish  at 
X  =  _+l/2.  Thus  'jj~ and  JJ”  are  determined  by  the  pair  of  differen¬ 
tial  equations 


subject  to 


TT 

+  <r  l,  + ■r- .  0 


rr,(±&)^Z(±'A)^o. 


(5.8) 


(5.9) 


(5.10) 


A  variety  of  elementary  methods  exist  for  solving  this  system. 
The  route  entailing  the  least  manipulation  probably  results  from  com¬ 
bining  (5.  8)  and  (5.  9)  into  one  second  order  equation  for  the  complex 
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variable 


'JJ^  -Xjf  Carrying  out  the  details  leads  us  to 


/  \  Jt  /  A  J 


(5.11) 


flfx) <Wi 4X dfio/  -  Sookt Oav i ,tutJu4  X tw4X 


(5.  12) 


where 


(5.  13) 


The  squeeze-film  force  is  given  by 

f‘/z 

F= J  CP-l)dX, 


(5.  14) 


where,  at  any  specific  time,  the  magnitude  of  F  is  the  average  loading 
pressure  (gage)  in  atmospheres.  Thus,  with  Eqs,  (5.  3)  and  (5.6), 

fivro  (5.  15) 

-'h 

With  (5.  11)  and  (5.  12), 


where 


FIs  -  (a-)  fiJoi  T+  ((r)xMx.  T 


icu^  ^/z 


(5.  16) 


(5.17) 
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^  (cr)  ~  jCU/n^ 


(UkK  -f  Qoq^  {s7z 


(5.  18) 


Equivalently, 

¥l{  =  -AU)c^[T+^iF)]  ,  l5->’) 

where 


(5.  20) 


(P  ~  ^i4<!i2nrv  {I'/j./fJ.  (5.21) 

For  small  values  ofcr-  ,  the  in-phase  and  out-of -phase  compon¬ 
ents  of  the  force  are  given  by 

6'^llZO  -h,.,  (5.22) 

<r^/lZ  -h  ,,,  .  (5.23) 

Thus,  the  amplitude  of  F /cT  for  smaller—  is  very  nearly  equal  to 
^2  {<’')  ih®  force  leads  the  motion  by  nearly  90°. 

For  very  large^j —  , 


4  (<3"}'-  f^/<r  \ 


so  that  the  force  and  the  motion  are  almost  in  phase. 

Figure  5.  1  illustrates  the  force-frequency  curve  in  the  range 
where  the  most  significant  behavior  is  observed;  the  solution  for  an 
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incompressible  film  is  included.  Figure  5.  2  illustrates  the  force - 
frequency  curve  over  the  full  range. 


VI,  SMALL,,  PERIODIC  VARIATION  OF  THE  GAP 
BETWEEN  PARALLEL  DISKS 

We  now  turn  our  attention  to  the  squeeze  film  between  two  paral¬ 
lel,  coaxial  disks  of  radius  B.  It  is  clearly  most  convenient  to  use  the 
polar  coordinate  system  introduced  in  Sec.  11.  The  projection  of  the 
bearing  on  the  X3  =  0  plane  is 

I  (6.1) 

Because  of  the  axial  symmetry,  the  pressure  between  the  plates 
is  governed  by  the  form  of  the  squeeze-film  equation  given  in  Eq.  (2.  10), 
subject  to  the  boundary  condition 

P(1,T)=  I,  (6.2) 

and  to  the  restriction 

P  (o,t)  finite .  (6-3) 

We  note  in  passing  that  the  pressure  is  governed  by  this  same 


equation  and  these  same  boundary  conditions  when  either  or  both  of 


Fig.  5.2.  Gas  Squeeze-Film  Force  Amplit-ude  and  Phase  for 
Infinifely  Long  Plane  Parallel  Plates,  O^o'5'lOOO 
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the  disks  rotate  about  their  mutual  axis.  This  can  be  seen  from  Eq. 
(2.8);  For  such  rotation,  steady  or  time-dependent,  and  Uj^ 
vanish.  Because  of  this  fact,  the  solution  to  Eqs.  (2.  10),  (6.  2)  and 
(6.  3)  for  small,  sinusoidal  variation  of  H  has  already  appeared  in  the 
literature.  At  the  Ninth  International  Congress  of  Applied  Mechanics 
(Brussels)  1956,  Professor  M.  Reiner  of  the  Israel  Institute  of  Tech¬ 
nology  demonstrated  an  instrument  for  measuring  the  non-Newtonian 
properties  of  air.  A  disk  6.  7  cm  ip  diameter  was  spun  at  7,  000  rpm 
opposite  a  fixed  stator  disk  0.  02  mm  away.  A  manometer  connected 
to  a  4-mm  hole  in  the  center  of  the  stator  recorded  pressures  which 
differed  qualitatively  from  the  predictions  of  Newtonian  flow  theory. 
Taylor  and  Saffman  [l]  investigated  the  possible  effects  of  engineering 
imperfections  in  Reiner's  instrument.  They  calculated  the  pressure 
perturbation  which  would  result  if  the  disks  were  not  quite  parallel, 
or  if  they  were  vibrating,  so  that  the  gap  between  them  varied  sinu¬ 
soidally  with  time. 

If  the  gap  between  the  disks  varies  according  to 

Ha  (f+  ^  (6-4) 

we  expect  that 

^  ,  (6.5) 

We  now  substitute  Eqs.  (6.4)  and  (6.  5)  into  Eq.  (2.  10),  neglecting 
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the  terms  of  the  second  degree  or  higher  in  and  obtain 

Taylor  and  Saffman  £7j  showed  that  a  periodic  solution  to  Eq.  (6.6), 
subject  to  Eqs.  (5.  2)  and  (5.  3),  is  provided  by 

(H]  )Oi<m  'T~j 


where 


and’ 


Ij]  ^Oh  (V^ 'ien>  Jc  ^&c  (/PR)  (6.9) 

The  net  force  W  acting  to  keep  the  disks  apart  is  given  by 

W  =  ;2Tr6^J  a  ~ pf^)0LR  >  {6.io) 


From  Eqs.  (6.  5)  and  (6.  7), 


(6.  11) 


w'=  % (rUR  -hzS J'  Rjf^  (r)ur. 


*  A  misprint  occurs  in  the  Taylor  and  Saffman  paper  at  this  point. 

The  second  of  their  equations  (15),  which  corresponds  to  our  Eq.  (6.  9), 
should  have  one  side  or  the  other  multiplied  by  -1. 
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To  evaluate  the  integrals  appearing  in  Eq.  (6.  11),,  we  make  use 
of  certain  identities  involving  Bessel  functions,  summarized  here  for 


convenience; 


I  ^  -X-)  j 

0'^ 

yc  {.i  , 


(6. 12) 


(6.13) 


(6.  14) 


^  I -i  X  .f-  ^  ^ 

o  -p 


(6. 15) 


(6. 16) 


(6.  17) 


(6.  18) 


Eqs.  (6.  17)  and  (6.  18),  which  are  the  ones  of  interest  to  us,  are 
derivable  from  Eqs.  (6. 12)  to  (6.  16). 


With  Eqs.  (6.  8)  and  (6.9),  Eq.  (6- 1 1)  becomes 


WA-  -  ^U)(iB^T  -i- 


<6.  19) 
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where 


^  / fl-)  ;  j  -  /?  ^  ^  /? '}-  ^  "^tA/ f  /^) 


!/?)*■+  (rfcLf^)'  (6.20) 

^  V?  ^ Vo)  f  '^Jse., 


Equivalently, 


{£rJ^[T+ ‘^(crjj  j 


(6. 21) 


(6.  22) 


where 


(6.23) 

(6.24) 


^  =■  OACtamj  i^^J^,)- 

The  behavior  of  g^  and  g^  for  small  values  of  (T" can  be  ascertained 
with  the  help  of  Eqs.  (6.  14)  and  (6.  16),  together  with  Maclaurin  series 

To(.^)=  +  F(iiF  “^ns/'r  ^(6-25) 


26) 


3/2 


By  setting  z  =  i  X.  <  obtain 
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■^^71  ■*■ 

Equations  (6.  20)  and  (6.  21)  then  give  us 

(f;=  ^Vw+,..  , 
q^O^)  '  (T  j?  +  •>.  ^  , 


(6.27) 


(6.  28) 


(6.29) 


(6.  30) 


(6.31) 

(6.32) 


To  investigate  the  behavior  of  gj^  and  for  large  values  of  0~'  , 
we  employ  the  asymptotic  formulas 

Ui7rz)^^dJsa^(2-/%7r/z-Tr/‘f)j  z- 


With  Eqs.  (6.  14)  and  (6.  16),  we  find,  after  some  lengthy  but 
straightforward  calculations. 


'Uu-p'^  iniri.  - 

(6.  34) 

-‘k  f*k 

'id  e  -tt/?), 

(6.  35) 

ioii  P  'V  iurtvLpit-  TTAA 

(6.  36) 

(6.37) 

Equations  (6.  20)  and  (6.  21)  then  yield,  for  largecT”'  , 
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l<f)  J  7^cr  ^  (6.38) 

✓V^  {zJ7  .  {6-39) 

It  is  not  surprising  that,  as  — -approaches  infinity,  so  tliat  the 
effect  of  leakage  becomes  less  and  less,  the  bearing  force  has  the  same 
functional  form  for  the  circular  disks  as  it  did  for  the  infinitely  long 
bearing.  That  the  coefficients  agree  is  coincidence:  li,  for  the  breadth 
B,  we  were  to  select  the  disk  diameter  instead  of  the  radius,  Eqs.  (6.  38) 
and  (6.  39)  would  undergo  an  obvious  modification. 
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